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Abstract 

We report a measurement of CP asymmetry parameters in the b ccd-transition-induced 
decays B^{B'^) — > J/i/; tt^. The analysis is based on a 140 fb~^ data sample accumulated at the 
T(45') resonance by the Belle detector at the KEKB asymmetric-energy e'^e~ collider. We fully 
reconstruct one neutral B meson in the J/ijj tt*^ final state. The accompanying B meson flavor is 
identified by its decay products. From the distribution of proper time intervals between the two B 
decays, we obtain the following CP violation parameters: 

Sj/^^o = -0.72 ± 0.42(stat) ± 0.08(syst) 
Aj/^^o = -0.01±0.29(stat)±0.07(syst). 

PACS numbers: 
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INTRODUCTION 



In the standard model (SM), the Kobayashi-Maskawa (KM) quark-mixing matrix[l| 
has an irreducible complex phase that gives rise to CP violation in weak interactions. 
In particular, the SM predicts large CP- violating asymmetries in the time- dependent 
rates of and decays into a common CP eigenstate /cp|_2|. In the decay chain 
T(4S') — > B^B^ fcpftag, where one of the B mesons decays at time tcp to a final 
state fcp and the other decays at time ttag to a final state Aag that distinguishes between 
B^ and 5°, the decay rate has a time dependence given by j3[ 

V{At) = ^ ^^'^^""' {l + 1 ■ [SfcP sin(AmdAt) + Af^p cos(AmdAt)] |, (1) 

where tbo is the lifetime, Am^ is the mass difference between the two B^ mass eigenstates. 
At = tcp — ttag, and the 6- flavor charge q = +1 (—1) when the tagging B meson is a 5° 
{B^). The CP- violating parameters Sf^p and Af^p are given by 

_ 2Xm(A) _ |AP-1 

^^"-|A|2 + 1' |A|2 + 1' 

where A is a complex parameter that depends on both the B^-B^ mixing and on the ampli- 
tudes for P° and B^ decay to fcp- To a good approximation in the SM, |A| is equal to the 
absolute value of the ratio of the B'^ —>■ fcp to B^ ^ fcp decay amplitudes. 

CP violation in neutral B meson decays involving the b ccs transition has been 
established through measurements of the CP-violation parameter sin20i by the Belle 0] 
and BaBar collaborations. Here, the SM predicts Sf^p = — ^/sin20i, where ^/ = +1(— 1) 
corresponds to CP-even (-odd) final states; and Af^^p = (or equivalently |A| = 1) for both 
b — i> CCS and the leading contributions to 6 — > ccd (e.g. the tree diagram). Hence, for 
fcp = J/i^ 71"°; which is a CP-even final state, Sj/^j^o becomes — sin20i if the tree diagram 
dominates. 

If penguin contributions or other contributions are substantial, a precision measurement 
of the time-dependent CP asymmetry in 6 — >■ ccd may reveal values for Sj/^t^o and Aj/^t^o 
that differ from what is expected. Measurements of CP asymmetries in 6 — > ccd transition- 
induced B decays such as B^ ^ J/ip tt^ thus play an important role in ascertaining whether 
or not the KM model provides a complete description of CP violation in B decays. 

A study of CP asymmetry in P° — > J/ip vr^ decays has been reported by the BaBar 
collaboration 7J. In this paper we report a measurement of time-dependent CP violating 
parameters in B^ —>■ J/ip 7r° decays using the higher statistics data accumulated by the Belle 
detector. 



DATA SAMPLE AND EVENT SELECTION 

The results presented here are based on a data sample of 140 fb~^ (corresponding to 
15.2x 10'' BB pairs) collected at the T(45') resonance with the Belle detector [1] at the KEKB 
asymmetric-energy e+e~ (3.5 on 8 GeV) collider |9j]. At KEKB, the T(4S') is produced with 
a Lorentz boost of jS'-f = 0.425 nearly along the electron beamline {z). Since the P° and P° 
mesons are nearly at rest in the T(4S') center-of-mass system (cms). At can be determined 
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from the displacement in z between the fcp and /tag decay vertices: At ^ (-^cp — -^tag) / (^ic = 
Az//37c. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a three- 
layer silicon vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of 
aerogel threshold Cherenkov counters (ACC), a barrel-like arrangement of time-of-fiight 
scintillation counters (TOF), and an electromagnetic calorimeter comprised of CsI(T/) crys- 
tals (ECL) located inside a super-conducting solenoid coil that provides a 1.5-T magnetic 
field. An iron fiux-return located outside of the coil is instrumented to detect mesons 
and to identify muons (KLM). The detector is described in detail elsewhere Q. 

Hadronic events are selected if they satisfy the following criteria: at least three recon- 
structed charged tracks; a total reconstructed ECL energy in the center of mass (cms) frame 
in the range between 0.1 and 0.8 times the total cms energy; an average ECL cluster energy 
below 1 GeV; at least one ECL shower in the region —0.7 < cos^^ < 0.9 in the laboratory 
frame; a total visible energy, which is the sum of charged track momenta and total ECL 
energy, exceeding 0.2 times the total cms energy; and a reconstructed primary vertex that 
is consistent with the known interaction point. After the imposition of these requirements, 
the efficiency for selecting B-meson pairs that include a. J/ip meson is estimated by Monte 
Carlo (MC) simulation to be 99%. To suppress continuum events, we require the event 
shape variable R2 to be less than 0.5, where R2 is the ratio of the second to the zeroth 
Fox- Wolfram moment (lol |. 

J /tp mesons are reconstructed via their decay into oppositely charged lepton pairs (e+e~ 
or fi'^fi~). Leptons are selected by starting with charged tracks satisfying \dz\ < 5 cm, where 
dz is the track's closest approach to the interaction point along the beam direction. For 
electron identification, the ratio between the charged track's momentum and the associated 
shower energy {E/p) is the most powerful discriminant. Other information including dE/dx, 
the distance between the ECL shower and the extrapolated track, and the shower shape 
are also used. Muons are identified by requiring an association between KLM hits and 
an extrapolated track. Both lepton tracks must be positively identified as such. In the 
e~^e~ mode, ECL clusters that are within 50 mrad of the track's initial momentum vector 
are included in the calculation of the invariant mass (Mee(7)), in order to include photons 
radiated from electrons/positrons. The invariant masses of e"'"e~(7) and /i^yU~ combinations 
are required to fall in the ranges —0.15 < (Mj/^ — Mee(7)) < +0.036 GeV/c^ and —0.06 < 
(Mj/^ — M^u ) < +0.036 GeV/c^, respectively. Here Mj/^ denotes the world average of the 
J/ip mass |l4| . 

Photon candidates are selected from clusters of up to 5x5 crystals in the ECL. Each 
photon candidate is required to have no associated charged track, and a cluster shape that 
is consistent with an electromagnetic shower. To select 7r° 77 decay candidates for the 
^ J/ip 7r° mode, the energy of each photon is required to exceed 50 MeV (100 MeV) 
in the ECL barrel (forward and backward endcap). vr^'s are formed from photon pairs that 
have an invariant mass in the range 0.118 to 0.15 GeV/c^. 

J/ip and vr*^ candidates are combined to select B candidates. The B candidate selection 
is carried out using two observables in the rest frame of the T(4S') (cms): the beam-energy 

constrained mass Mbc = i/-Sbeam ~ (SPi)^ the energy difference AE = J^Ei — -Ebeam, 
where -Ebeam = is the cms beam energy, and pi and Ei are the cms three-momenta and 

energies of the B meson decay products. In this calculation, kinematic fits are performed 
with (1) vertex and mass constraints for the J/ip di-lepton decays and (2) mass constraint for 
the tt'^ — i> 77 decays in order to improve the AE resolution. Candidate events are selected by 
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requiring 5.270 < Mbc < 5.290 GeV/c^ and -0.10 < < 0.05 GeV. The lower bound of 
the /S.E requirement is determined in order to accommodate the negative /S.E tail that results 
from shower leakage associated with the high-momentum tt^. The number of reconstructed 
— > J/ip 7r° candidates is 103. The /S.E and Mbc distributions for the candidate events 
are shown in Fig, ^ 
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FIG. 1: The Aii^(left) and Mbc(right) distributions for J vr*^ candidates. The superimposed 

curves show fitted contributions from signal (red two-dot-dash, plotted only in Ai? distribution), 
B J/il:X background (green dot-dash), combinatorial background (blue dash) and the sum of 
all the contributions (black solid). See text for further details. 



SIGNAL PROBABILITY 

To assign an event-by-event signal probability for use in the maximum-likelihood fit to the 
CP-violating parameters, we determine event distribution functions in the Aii^-Mbc plane 
for both signal and background. The signal distribution is modeled with a two-dimensional 
function which is Gaussian in Mbc and uses a Crystal Ball line shape |13] in AE. The shape 
parameters of these functions are determined from MC simulation and held fixed in the 
fit, while the overall signal yield is allowed to float. Backgrounds are studied using a large 
sample of MC events along with events outside of the signal region. We split the backgrounds 
into two categories, one being B decays having a J/ip{B J/ipX) and the other being 
combinatorial background to which random combinations of particles in BB decays and 
continuum events contribute. According to MC study, the B —>■ J/ipX background forms a 
small peak in the Mbc projection. Therefore we parametrize this contribution with the sum 
of a Gaussian and a phase-space like background function(ARGUS function) |1^ in the Mbc 
direction and an exponential function for AE. The amount of this background contribution 
is determined by MC |17|]. For the combinatorial background, we use a linear function for 
AE and an ARGUS function for Mbc. The purity of the signal is estimated to be 86±10% 
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FLAVOR TAGGING AND VERTEXING 



Charged leptons, kaons, pions, and A baryons that are not associated with the recon- 
structed 5° — >• J/V' vr" decay are used to identify the 6-flavor of the accompanying B meson, 
denoted by /tag- Based on the measured properties of these tracks, two parameters, q and r, 
are assigned to each event. The first, q, has the discrete value +1 (—1) when the tag-side B 
meson is more hkely to be a B^ (-B")- The parameter r is an event-by-event MC-determined 
flavor-tagging dilution factor that ranges from r = for no flavor discrimination to r = 1 
for an unambiguous flavor assignment. It is used only to sort data into six intervals of 
r, according to the estimated flavor purity. The wrong-tag probabilities for each of these 
intervals,w/ (/ = 1,6), which are used in the final fit, are determined directly from the data 
samples of B^ decays to exclusively reconstructed self-tagging channels. The difference of 
the wrong-tag fractions between B^ and B^ are also determined from the same samples and 
implemented as Awi{l = 1, 6) into the final fit. We obtain wi and Awi using time-dependent 
mixing: (A^of ~ -^sf)/(-^of + -^sf) = (1 2w/) cos(AmAt), where A^of and Ngp are 
the numbers of opposite B^B^) and same B^B°, B^B^) flavor events. 

The wrong tag fractions for each r interval are given elsewhere jl^ . 

The decay vertices of B^ mesons are reconstructed using tracks that have enough SVD 
hits: i.e. both z and r-0 hits in at least one SVD layer and at least one additional layer 
with a z hit, where the r-0 plane is perpendicular to the z (beams) axis. Each vertex 
position is required to be consistent with the IP profile, which is determined run-by-run and 
smeared in the r-0 plane by 21 /xm to account for the B meson decay length. With these 
requirements, we are able to determine a vertex even in the case where only one track has 
enough associated SVD hits. The vertex position for the B^ J/tp ir^ decay is reconstructed 
using lepton tracks from the J/ip. The algorithm for the /tag vertex reconstruction is chosen 
to minimize the effect of long-lived particles, secondary vertices from charmed hadrons and a 
small fraction of poorly reconstructed tracks . From all the charged tracks with associated 
SVD hits except those used for B^ ^ J/ip 7r° reconstruction, we select tracks with a position 
error in the z direction of less than 500 /im, and with an impact parameter with respect 
to the J/ip vertex of less than 500 fim. Track pairs with opposite charges are removed if 
they have an invariant mass within ±15 MeV/c^ of the nominal Ks mass. If the reduced 
associated with the /tag vertex exceeds 20, the track making the largest contribution is 
removed and the vertex is refitted. This procedure is repeated until an acceptable reduced 

is obtained. After flavor tagging and vertex reconstruction, we obtain 91 B^ ^ J/ip 7r° 
candidates. 



THE UNBINNED MAXIMUM LIKELIHOOD FIT 

We determine Sj/^t^o and Aj/^-j^o for each mode by performing an unbinned maximum- 
likelihood fit to the observed At distribution. The probability density function (PDF) 
expected for the signal distribution is given by 

Vsig{At, q,wi,Awi) 

= - + - 2^') ■ ['^/./,.o sin(AmrfAt) + Af^^^^^, cos(Am,At)] | (3) 

to account for the effect of incorrect flavor assignment. The distribution is convolved with the 
proper-time interval resolution function Rsig{At), which takes into account the finite vertex 
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resolution. i?sig(At) is formed by convolving four components: the detector resolutions for 
zcp and ^tag, the shift in the Ztag vertex position due to secondary tracks originating from 
charmed particle decays, and the kinematic approximation that the B mesons are at rest 
in the cms A small component of broad outliers in the Az distribution, caused by 

mis- reconstruction, is represented by a Gaussian function Poi{At). We determine twelve 
resolution parameters and the neutral- and charged-5 lifetimes simultaneously from a fit to 
the At distributions of hadronic B decays and obtain an average At resolution of ~ 1.43 ps 
(rms). We determine the following likelihood value for each event: 

/oo r 
/sigPsig(At', q, wu Awi)Rsig{Ati - At') 
-oo L 

+ (1 - /sig - fiir)nifWK^fiAt, - At')] d{At') + /oiPoi(At.) (4) 

where foi is the outlier fraction and /sig is the signal probability calculated as a function of 

AE and M^c- vH^^ (At) and V^lf^{At) are the PDFs for B J/ipX and combinatorial 
background events, respectively. These contributions dilute the significance oi CP violation 
in Eq. (^. They are modeled as a sum of exponential and prompt components, and are 
convolved with the corresponding resolution functions -Rb^kg^ ^^"^ Phkg^^ respectively. The 
resolution functions are modeled by a sum of two Gaussians. All parameters in V^lg^ {At) 
and -Rbkg^ determined from MC simulation, while the parameters in V^^'^^{At) and 



R^^^ are determined by a fit to the At distribution of a background-enhanced control 



■-bkg 

■-bkg 

sample; i.e. events away from the AE-M]^c signal region. We fix r^o and Arrid at their 
world-average values [1J|. The only free parameters in the final fit are Sj/^^^o and v4j/^^o, 
which are determined by maximizing the likelihood function 



^ = '[[PiiAti; Sj/^^o , Aj/^T,o ) (5) 

i 

where the product is over all events. 



FIT RESULTS 

A fit to the candidate events results in the CP-violation parameters; 

Sj/^no = -0.72 ± 0.42(stat) ± 0.08(syst) 

Aj/^^o = -0.01 ±0.29(stat) ±0.07(syst), (6) 

where the sources of systematic error are given below. Figure |2] shows the At distributions for 
B^ ^ J/ip 7r° (upper figure:g = -|-1) and B^ ^ J/ip 7r° (lower figure:g = —1) event samples. 
Figure 01 shows the raw asymmetry in each At bin without background subtraction, which 
is defined by 
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At(ps) 

FIG. 2: The At distributions for B° J/ip 7r° (upper:g = +1) and Jj^p 7r° (lower :g = -1) 

candidates. The sohd curves show the results of the global fits, and dashed curves show the 
background distributions. 

where N q=j^\{N q=-\) is the number of observed candidates with q = 1). The 

curve shows the result of unbinned-maximum hkehhood fit to the At distribution, 
Sj/^^o sin(AmdAt) + Aj/^^o cos(AmdAt). 

SYSTEMATIC UNCERTAINTIES 

We estimate systematic uncertainties as follows; the fiavor tagging (±0.025 for Sj/^-^o 
and ±0.014 for ^j/^^o), the signal probability (±0.023 for Sj/^^^o and ±0.016 for Aj/^tto), 
the background At distribution (±0.014 for Sj/^^^o and ±0.0065 for Aj/^^^o), the resolution 
function (±0.010 for Sj/^t,o and ±0.0061 for Aj/^t^o), the potential fit bias (±0.043 for 
Sj/^ttO and ±0.059 for Aj/^t,o), the vertex reconstruction (±0.062 for iSj/^^o and ±0.018 for 
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FIG. 3: The Ai asymmetry. The curve shows the result of the unbinned-maximum hkehhood fit. 



Aj/^tto), and B meson's lifetime and mixing parameter (±0.0019 for Sj/^^^o and ±0.0071 for 
Aj/^tto)- The quadratic sum of all the contribution mentioned above amounts ±0.0845 for 
Sj/^-nO and ±0.0662 for Aj /^t^o. 



SUMMARY AND CONCLUSIONS 

We have performed a measurement of CP- violation parameters in —>■ J/ip tt^ de- 
cay. The resultant values are iSj/^^o = —0.72 ± 0.42(stat) ± 0.08(syst) and Aj/^ttO = 
—0.01 ± 0.29(stat) ± 0.07(syst). These values are consistent with those obtained for 
J/ip Ks and other decays governed hy h ccs transition and suggest that pen- 
guin and other contributions to this decay mode are not large. 
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